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Doublon-hole pair production which takes place during dielectric breakdown in a Mott insulator 
subject to a strong laser or a static electric field is studied in the one-dimensional Hubbard model. 
Two nonlinear effects cause the excitation, i.e., multi-photon absorption and quantum tunneling. 
Keldysh crossover between the two mechanisms occurs as the field strength and photon energy 
is changed. The calculation is done analytically by the Landau-Dykhne method in combination 
with the Bethe ansatz solution and the results are compared with those of the time dependent 
density matrix renormalization group. Using this method, we calculate distribution function of the 
generated doublon-hole pairs and show that it drastically changes as we cross the Keldysh crossover 
line. After calculating the tunneling threshold for several representative one-dimensional Mott 
insulators, possible experimental tests of the theory is proposed such as angle resolved photoemission 
spectroscopy of the upper Hubbard band in the quantum tunneling regime. We also discuss the relation 
of the present theory with a many-body extension of electron-positron pair production in nonlinear 
quantum electrodynamics known as the Schwinger mechanism. 
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I. INTRODUCTION 

"Nonequilibrium strongly correlated systems" is be- 
coming an important field of study in condensed matter 
physics p]4l^|. These systems offer a tcstbed for theo- 
retical advances such as the extension of the linear re- 
sponse paradigm to nonlinear processes fl3rl26l ]. We can 
experimentally induce a nonequilibrium state in photo- 
induced phase transitions in solids 0, Q as well as in 
the dynamics of cold atoms [9rJl2j]. The photo- induced 
insulator to metal transition in Mott insulators has gen- 
erated substantial interest because it is one of the most 
basic nonequilibrium phenomena in strongly correlated 
systems [l|, [|[ . The response of Mott insulators sub- 
ject to strong external fields has been studied experi- 
mentally. Initially, the motion of particles (electrons or 
atoms) is frozen by strong repulsion, and the ground state 

Perturbations (electric field or 



extensively in theory via the fermionic Hubbard model 
usin g n umerical methods such as exact diagonalization 
fl3l |21| , the time-dependent density matrix renormal- 
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is a Mott insulator 
lattice modulation) excite pairs of doublons (= doubly 
occupied site) and holes (= sites with no electrons; we do 
not call this state "holon" because this is a Bethe ansatz 
terminology that is used later), and when their density 
becomes sufficiently high, "melting" of the Mott state oc- 
curs leading to an insulator-to- metal transition [l|. Quite 
recently, the insulator-to-metal transition was realized by 
a terahertz laser in vanadium dioxide, which is a candi- 
date material for a Mott transition Q . Since the photon 
energy is far below the optical gap, the excitation mech- 
anism is expected to be a nonlinear process. These ex- 
periments give us strong motivation to develope a theory 
for nonlinear excitations in strongly correlated systems. 

The purpose of this study is to gain an analytical un- 
derstanding of the excitation process when a strong elec- 
tric field is applied to a Mott insulator. The effects of 
strong electric fields on Mott insulators have been studied 



2C|, and nonequilib- 
17, [H. These stud- 



ization group (td-DMRG) 0, H 
rium dynamical mean field theory 
ies reveal the following consensus. Doublon-hole pairs 
(dh-pairs) are created by strong electric fields, and for 
DC-electric fields, production rate (or ground state de- 
cay rate) shows a threshold behavior [t| [l4|. This be- 
havior seems to be universal and independent of dimen- 
sions, e.g., Ref s. fl3Ml5l [20l| (one-dimensional (1-D) stud- 
ies) and Ref. (infinite-dimensional studies). If we 
denote tunneling threshold by F t h, for small electron re- 
pulsion U, it behaves as F t h oc A Mott , where A Mott is 
the Mott gap. We can obtain this expression by ap- 
plyin g th e Landau-Zencr formula to many-body energy 
levels 1131 . For AC-clectric fields, it was mentioned in 
Ref. [2l[ that there is a crossover from a weakly ex- 
cited state to a strongly excited state with increasing 
field strength. Another interesting observation was made 
regarding the bosonic Hubbard model with lattice mod- 
ulation, where the authors calculated energy absorption 
rate using td-DMRG [l2[ . The absorption peaked around 
51 ~ NU (CI: modulation frequency, N: integer) broad- 
ened as modulation intensity increased. The broadening 
is clearly a nonlinear effect. 

In this study, we examine the ID Hubbard model 
at half-filling and the instability of the ground state in 
strong electric fields. The method we use is a combi- 
nation of the Landau-Dykhne quantum tunneling the- 
ory (28rl30| and the Bethe ansatz. This method was 
developed in Ref. [l5j and was used to derive an ana- 
lytic expression for tunneling threshold. Although not 
commonly employed in condensed matter, the Landau- 
Dykhne method has a long application history in areas 
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of physics where quantum systems are driven out of their 
initial state by strong external fields. The name "strong 
field physics" is often used to describe this problem area 
in physics. We can find examples of driven systems in 
quantum chemistry (30j , atom ionization (3ll [32j , quan- 
tum chaos [33^ . and high energy [33 - l4ll |. "Nonlinear 
excitations in Mott insulators" is a typical problem in 
"str ong field physics in condensed matter" (for a review, 
see [U HU). Ideas and techniques developed in other 
fields prove quite useful as well. 

Our problem has many common features with the 
electron-positron pair production problem in nonlin- 
ear quan tum electrodynamics (QED) (for a review, see 
Ref . [4l[ ) . The concept of vacuum in high energy physics 
is directly translated into an "insulator" in condensed 
matter. Shortly after the report of the Dirac sea vac- 
uum description, Heiscnbcrg-Euler proposed that non- 
linear response of the vacuum is described by an ef- 
fective Lagrangian (34|. They also found vacuum in- 
stability against electron-positron pair production when 
field strength is comparable to tunneling threshold. This 
threshold is now called the Schwingcr limit [35|. The 
calculation of production rate was extended from DC 
to AC electric fields [H-HH following an early study by 
Keldysh on atom ionization j3l| . In particular, Popov 
used the Landau-Dykhne method to calculate production 
rate [36|, H3] (for a review of this approach, see Ref. [39| ) . 
Following these "strong field physics" -studies, a universal 
picture emerged that, in fact, had already been noticed 



by Keldysh[31|. That is, there are two leading excitation 
mechanisms in a zero-temperature-gapped system driven 
by an AC-ficld. One is quantum tunneling, which is dom- 
inant in the DC limit, and shows a threshold behavior. 
This threshold is nothing but the Schwinger limit. The 
other mechanism is multi-photon absorption that is dom- 
inant when photon energy is relatively large. Moreover, 
production rate shows a power law behavior. There is 
a crossover between the two regimes, which is called the 
Keldysh crossover. 

In this study, we show that the nonlinear dh-excitation 
in a Mott insulator also lies within the Keldysh paradigm. 
Using the Landau-Dykhne method combined with the 
Bethe ansatz, we derive an expression for the momentum- 
resolved production rate of dh-pairs (Eq. below) 
and calculate the total production rate T. In Fig. [TJ we 
schematically plot the production rate behavior in strong 
AC-clcctric fields 



F(t) = F sin nt. 



(1) 



Here, Fq = eaE(> 0) is the field strength, e is the elec- 
tron charge, a is the lattice constant, and f2 is the photon 
energy. When photon energy is above the Mott gap and 
is resonant with the absorption spectrum, i.e., fl ~ U, we 
obtain the standard linear response result, i.e., T oc Fq. 
In the case where the field is off-resonant 51 < Amou, 
nonlinear processes lead to dh-production. Similar to the 
other above-mentioned "strong field physics" examples, 
the two leading mechanisms arc multi-photon absorption 
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FIG. 1: (color online) (a) Schematic plot of the nonlinear 
optical absorption spectrum of a Mott insulator in a strong 
AC-electric field. Fo(> 0) is the field strength and f2 is the 
photon energy of the applied laser. In addition to the contri- 
bution from linear response theory (Kubo formula), sub-gap 
excitations occur owing to nonlinear processes. Mechanisms 
are quantum tunneling and multi-photon absorption. They 
are governed by the Mott gap Aiviott and correlation length 
£. (b) Correlation length f is the typical size of doublon-hole 
pairs in the Mott insulating ground state. Pair production 
and annihilation occur during the virtual process. 



and quantum tunneling. The former occurs when is 
relatively close to the gap. Production rate has the fol- 
lowing power law dependence on field strength 



roc (H) 



(2) 



Here, the power 2 



A 3 



is twice the number of absorbed 



photons and the factor £ is the doublon-hole correlation 
length [45| . In the ground state of a Mott insulator, 
doublon-hole pairs are created during a quantum me- 
chanical virtual process [Fig. [Hb)]. Correlation length 
gives the typical size of such doublon hole pairs. When 
the DC limit is approached with a small f2, the leading 
mechanism becomes quantum tunneling. This leads to 
a dielectric breakdown with a threshold behavior (]~3l — 
EE [13 ■ The total production rate in this regime has the 
approximate form 



r oc exp ( — 7T— - 



F 



Hi 



(3) 
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where DC threshold (= Schwingcr limit) is given by 



by[4J 



A 



Molt 



2£ 



(4) 



We notice that the correlation length £ again plays an 
important role. 

The paper is organized as follows. In section HH after 
a brief introduction of the Bethe ansatz solution of the 
Hubbard model, we explain how to combine its infor- 
mation with the Landau-Dykhnc method. Application 
to nonlinear transport in DC fields and photo-induced 
phase transitions in AC fields is discussed in sections IIIII 
and llVl respectively. In section [Vj we discuss experimen- 
tal feasibility. 



II. LANDAU-DYKHNE + BETHE ANSATZ 
METHOD 



In this section, we extend the Landau-Dykhnc + Bethe 
ansatz method, developed in Ref. [l5j], to electric fields 
with various laser types. The model we study is the half- 
filled ID fcrmionic Hubbard model subject to an electric 
field. The Hamiltonian is given by 

H(*) = -ryej^ + e^i^) (5) 

3 

The time-dependent Peierls phase <& is related to the ap- 
plied electric field by F(t) = eaE(t) = -d$(t)/dt. We 
set the energy unit as the hopping amplitude, i.e., r = 1. 
We start from the Mott insulating ground state at t = 
and apply the electric field for t > 0. 

The static Hubbard model can be solved exactly us- 
ing the Bethe ansatz and the ground state wave func- 
tion as well as excitations is well understood [44| . There 
arc two types of elementary excitations from the half- 
filled ground state: (1) Gapped spinlcss excitations 
with charge =Fe called antiholons and holons (2) Gap- 
less charge neutral excitations carrying spin ±i called 
spinons. Physical excitations are built from these ele- 
mentary excitations. In the remainder of this article, 
instead of using the Bethe ansatz terminology antiholon 
and holon, we use the more familiar names doublon and 
hole. 

Among the excitations, we concentrate on excited 
states with a single doublon- hole pair, i.e., antiholon- 
holon pairs. The states are parameterized by rapidity k\ 
(fo) for a hole (doublon). The total energy and central 
momentum of these excitations are given by 

AE = E dh -E = e h {kx) + e d {k 2 ), (6) 

-Pccntral = Ph(fa) + Pd(fe), (7) 

where Eh,d and ph,d are the energy and momentum of the 
hole and doublon, respectively. Holon energy is given 



e d (k) = 172 + 2 cos fc 

doj Ji(w) cos(w sin k)e~ UuJ / 4 



+2 



OJ 



cosh(w[//4) 



(8) 



and holon momentum is given by 



Ph(k) = Pd{k) + n = — - k 



doj Jq (uj) sin(o; sin(fc) ) 
bj 1 + cxp(J7o;/2) 



(9) 



Note that we shifted holon momentum by ^, i.e., p = 
Ph + f (= ~Pd - § )• Only states with P CG ntrai = can 
be excited by an external electric field because the mo- 
mentum of the laser can be ignored. We denote a single 
doublon- hole pair with the hole momentum p by \p)dh- 
In Fig. HJa), we plot the excitation energy AE(p) as a 
function of p. It has a minimum at p = with a gap 
A Mo tt, i.e., the Mott gap. 

An important concept in the Mott insulating phase is 
the correlation length £ studied by Stafford and Millis 
in Ref. [45| using the Bethe ansatz. In Mott insulators, 
each lattice site is occupied by a single electron. How- 
ever, quantum fluctuations enable doublons and holes to 
pair create and wander around as a virtual process be- 
fore they pair annihilate. This process is responsible for 
antifcrromagnetic super exchange coupling. Correlation 
length is intuitively the size of the doublon-hole pair in 
the ground state wave function [Fig.fljb)]. If we consider 
a finite system of size L, charge stiffness (Drude weight) 
decays as 



D C (L) oc exp(-L/£). 



(10) 



In other words, if the system is small enough compared 
to £, it behaves as a metal because the carriers (dou- 
blons and holes) can transport current. Another im- 
portant property is that Green's function G(\x — x'\) = 



decays as follows [4£ 



G(|x|)~cxp(-M/0- 
The exact expression for £ is given by (45 



4 p ln(y +v ^i) 
/? U J 1 V cosh(2Try/U) ' 



(11) 



(12) 



In Fig. EJa), we plot £ as well as Amou as a function of 
U. In the small U limit, it behaves as 



lim £ = 



2t + U/2-k ■ 



(13) 



A(U,t) 

whereas the strong-coupling limit is given by 

C 1 = HU/at) U — y oo (14) 
with a = [r(l/4)/ % /2^] 4 ~ 4.377. 
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FIG. 2: (color online) (a) The excitation energy AE of the doublon-hole pair as a function of the hole momentum p in the 
U — 8 Hubbard model. V v is the tunneling probability to create a state \p)dh from |0). (b) The real part of the excitation 
energy plotted for complex p. The Imp = slice is equivalent to the left half of (a) . The gap closes at the level crossing point 
p — ip c where a gapless line starts. Paths Fi and F2 are used in the integral in Eq. {25) . 




FIG. 3: (color online) (a) Correlation length £ [45| and (b) 
Mott gap Aiviott (Lieb-Wu solution) of the ID Hubbard model 
at half-filling. 



Next, we consider time evolution. We start from the 
ground state |0) and apply an electric field described by 
a time-dependent phase via F(t) = eaE(t) = —d§(t)/dt. 
After a tunneling process, the wave function takes the 
form 



I*) 



)dh + • ■ 



(15) 



(a, ftp are phases) where the ground state amplitude de- 
creases as VJ7 • 



lip \/l — ■ 'Pp is the momentum 



with momentum p [Fig. [2] (a)]. The omitted term ". . ." 
contains excitations to states with multiple doublon-hole 
pairs as well as spin excitations. 

One can calculate the tunneling probability P n in the 
Hubbard model by the Landau-Dykhne method [15[ . The 
Landau-Dykhne method [2£jL l3(i| (for a textbook and use- 
ful reference see Ref. [32| and j46l| , respectively) has been 
derived from the adiabatic perturbation theory. We de- 
note the adiabatic eigenstates of H(&) by |0;$) and 
\p;®)dh where 

£T($)|0;$) = £ ($)|0;$>, (16) 
H($)\p;$) dh = E dh (p;$)\p;$) dh (17) 

is satisfied. Because p is a good quantum number, states 
with different p are orthogonal to each other. Thus, by 
ignoring multiple pair states, we can study each exci- 
tation channel independently. This means that we can 
study the tunneling process in a Hilbert space spanned 
by two states, i.e., the problem reduces to solving the 
time-dependent Schrodinger equation with a solution of 
the form 



\*(t)) = a(t)\0;$(t)) +b(t)\p;$(t)) dh 



(18) 



(initial condition a(0) = 1, 6(0) = 0). This significantly 
simplifies the problem. Landau-Dykhne's tunneling the- 
ory states that the tunneling probability (|6(i C nd)| 2 ) be- 
tween two quantum levels is given by 



where 



7> p = exp(-2Im2? p ), 



V P = [ [E dh (p;<f>(t))-E ($(t))]dt 



(19) 



(20) 



resolved tunneling probability of the doublon-hole pair 



is the difference between the dynamical phase of the 
ground state and the excited state. In our problem, 
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this expression can be simplified because the effect of 
the Peierls phase on the adiabatic solutions is expressed 
simply by replacing the momentum p by p— $. Thus, we 
have 



llh 



\p - $} dh , E dh {p- $) = E dh (p - (21) 



which leads to V p = f AE(p - $(t))dt, where AE is 
defined in Eq. flBJ). 

An interesting point of the Landau-Dykhnc formula 
[Eq. (|19p] is that tunneling probability depends on the 
imaginary part of the dynamical phase difference. Inte- 
gration path 7 starts from t = and ends at a critical 
time t = t c , where level crossing takes place. Because 
we are dealing with a gapped system, level crossing does 
not occur for real t; instead, it occurs at a complex time 
when 



AE(p - $(t c )) = 



(22) 



is satisfied. When t is a complex number, the correspond- 
ing Peierls phase is also complex. The Hubbard Hamil- 
tonian [Eq. ((SJ)] with a complex Peierls phase is a non- 
Hermitian lattice model where the absolute values of the 
left and right hoppings are unequal. The ground state 
wave function of the non-Hcrmitian model was studied 
by Fukui and Kawakami in ref. [47| . In Fig. GJb), we 
plot TLeAE for complex p. The level crossing is found 
at a point ip c = p — &{t c ). The momentum of the level 
crossing point is related to correlation length by [Ifl [48{ 



pc = ve, 



(23) 



i.e., they are the inverse of each other. In a nonintcr- 
acting system, this is a very natural relationship, which 
states that localization length is the inverse of complex 
momentum and a wave function decays as e lkx = e~ KX 
when k = in. However, the surprise here is that this 
concept can be extended to a many-body system in a 
straightforward manner. 

The expression for tunneling probability becomes 
physically clearer when we change variables in the in- 
tegral from time t to the Peierls phase Using the 
Jacobian ^ = —1/F, where F is the electric field ex- 
pressed as a function of <f>, we are led to the expression 



cxp -21m / AE(p - <f>) 



■1 



F($) 



(24) 



which is the main result of this study. We keep the mi- 
nus sign in Eq. (|24[) as a reminder that the factor in the 
exponential is negative. This expression is a direct de- 
scendant of V. S. Popov's original expression for the tun- 
neling probability in the massive Dirac model[37|. The 
difference here is that the one-body energy level is re- 
placed by the many-body level obtained by the Bethe 
ansatz. One can use Eq. (|M| to study excitations not 
only for DC-fields but also for various other fields. This 
can be achieved by simply replacing the function F(<f>) 



TABLE I: Models of electric fields. F is the field strength, fi 
the photon energy and a the pulse duration. 



type 


F(t) 


F($) 


attempt frequency / 


DC-field 


F 


F 


Fq/2-k 


AC-field 


Fo sin Qt 


±V F o - ft 2 $ 2 


Q/2tv 


single pulse 


F cosh -2 (t/a) 


*('-*) 


1 (single process) 



in the formula. Table U shows the models of the electric 
fields we use in this study. 

There is an arbitrariness in the $ integral in Eq. (|24l) . 
Because T> p is a complex integral of an analytic function, 
one can deform its path as long as the end points are 
fixed and no singular point is crossed. The most natural 
path that simplifies the calculation is the one that goes 
from p to (ri) and then from to ip c (T2) as shown in 
Fig- [2Kb). On this path, AE is always real. The integral 
is divided into two 



<iv 



ae{ p - $; 



d$ = V, 



V 



p2, 



(25) 



corresponding to Y\ and T2, and the two contributions 
can be written as real integrals 



InrP, 



pi 



AE(l)lm 



1 



F(p - I) 

LmP p2 = /q AE(U)Jm 



dl, (26) 
dl. (27) 



In these integrals, the variable I is real. We can numeri- 
cally perform the integration, and in some simple cases, 
we can derive analytical expressions, as we will see later. 

We now comment on the relationship between tunnel- 
ing probability and production rate. The tunneling prob- 
ability V p is defined for a single excitation attempt. To 
promote it to a "rate" , we must multiply it with the 
attempt frequency /, i.e., the number of events per unit 
time. Then, momentum resolved production rate is given 
by 



(28) 



In the case of DC fields, physical momentum, defined on 
the interval (— tt, it] , evolves as p — <3?(i) = p + Fat. Thus, 
the time period between tunneling events is 2tt/Fq:, 
which is nothing but the period of the Bloch oscillation. 
In this case, the attempt frequency is given by its inverse 
/ = i<o/27r. As for the AC field with photon energy ft, 
the attempt frequency is / = 0/2vr. When we consider 
a single pulse, we simply set / = 1. Table U summarizes 
attempt frequency and the Jacobian for several models 
of electric fields. 

The total production rate is defined by 



r-fl 



(29) 



The total production rate is an important quantity be- 
cause it is comparable to quantities obtained from with 
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other methods such as td-DMRG. First, the total pro- 
duction rate gives the lowest order approximation for 
the ground state decay rate 0]. The ground-state- 
to-ground-state transition amplitude (fidelity amplitude) 
for a time-dependent Hamiltonian is defined as 



B(t) = (0;$(T)|Te-^° ff( * (s))ds |0;$(0)) 



(30) 



where T stands for time ordering. When the ground state 
is unstable in the external driving force, the absolute 
value of the amplitude decays exponentially as (D = 1 
ID) 



|S(f)| 



-tL u V„ 



(31) 



The total production rate defined in Eq. (|3"Tj) agrees with 
the ground state decay rate up to the higher order tun- 
neling process, i.e.,r g s . decay = T + . . ., where terms such 
as (Tp) 2 are neglected. For example, the ground state 
decay rate of a band insulator in DC fields is given by 

r g . s .dccay = -*7 B .z. (JrF ln i 1 ~ V k) WnC1 ' C fc is tllC 

momentum in the Brillouin zoneflij. Expanding this 
with V to the lowest order gives the total production 
rate of electron- hole pairs, c.f. Eq. (|2T))) . 
The time evolution of the doublon density 



1 L 

d(t) = - 2j(n it n i4 .)(t). 



(32) 



can be related to the total production rate. As we apply a 
strong electric field to the ground state, doublon density 
increases from its ground state value. For a continuously 
applied electric field, we assume that doublon density 
increases linearly in time, and its increase within a time 
interval At is given by 



Ad ~ Air. 



(33) 



Again, this expression is approximate because (a) we ig- 
nore the production of multiple pairs and (b) we assume 
that the state \p; $) has precisely one additional doublon 
compared with the ground state, which is a natural as- 
sumption when U is large. 

Next, we comment on the validity of our method. It is 
important to note that although we use an exact result 
(the Bethe ansatz), the calculated production rate is only 
approximate. One origin of error lies in the Landau- 
Dykhne formula itself. It is only valid when excitations 
are rare events, i.e., V p <C 1. This means that we can only 
use Eq. ([24)) when field strength is not too large compared 
with the Schwinger limit F t h and photon frequency is 
below the resonance frequency f2 < Amou- A related 
issue is that the Fq ~ U resonance [24| is ignored in 
static electric fields. 

Another important omission is the effect of quantum 
interference between multiple tunneling events. This is 
known as the Stokes phenomenon and has been studied in 



various time-dependent problems (e.g., Ref. (la. liol. l46j ) . 
In the present problem of a driven Mott insulator, this is 
related to the pair annihilation process of doublon- hole 
pairs. In Ref. [161, the effect of 

pair annihilation and the 
resulting quantum interference was studied via mapping 
into an effective quantum walk. Quantum interference 
may lead to several anomalous behaviors. An example is 
dynamical localization in energy space [l6| . The Landau- 
Dykhne method ignores the interference effect, and we 
will see the outcome of this later in Section IIVBI while 
presenting a comparison with numerical results. 



III. DIELECTRIC BREAKDOWN IN DC 
FIELDS (SCHWINGER LIMIT) 



"Upper Hubbard band" 




'jwer Hubbard band" 



A 



Mott 



FIG. 4: (color online) Schematic rigid band description of the 
dielectric breakdown of Mott insulators in DC-electric fields. 
The upper and lower "Hubbard bands" are tilted with static 
potential V(x) = Fqx and quantum tunneling starts to occur 
when the energy drop Fo(, between the doublon pairs sepa- 
rated by £ becomes comparable to the excitation gap Amou • 

The case of DC electric fields was studied in Ref. [l5| , 
and the Schwinger limit (=tunncling threshold) was 
obtained and compared with td-DMRG. In DC fields 
F = Fq, when we perform the integrals in Eqs (|26[) 
and (|27|) . we notice that the contribution from path 
Ti vanishes because F is always real. Thus, we have 
ImVp = Jq° AE(il)dl/Fo which leads to the threshold 
form 



exp 



.Fth 
F 



with 



F 



th 



Pc=l/£ 



AE(il)dl. 



(34) 



(35) 



This coincides with earlier results [15J. Further physical 
insight can be gained using an approximation (accurate 
for small U I50l0 such as 



AE(il) ~ A M ott\/r 

which makes the integral in Eq. 
obtain the Schwinger limit 

Amou 



■ (CO 2 , 

35l) trivial. 



F th 



2£ 



(36) 
Then, we 

(37) 
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As explained above (section [TTJ) , attempt frequency is 
given by 1/T = Fq/2t:, and production rate becomes 



r F o ( -Fth 
r^-expl-vr— 



(DC-fields). 



(38) 



We note that the case of a DC field is special in the 
sense that production rate has no p-dcpcndcncc. This is 
because all states experience the same tunneling event. 
A state with momentum p drifts in the momentum space 
as p + F$t and undergoes tunneling around p + F$t ~ 
when the gap becomes smallest. 

The interpretation of the result can be simplified if we 
employ the rigid band picture (Fig. H|) . The rigid band 
picture simply views a Mott insulator as a band insulator 
with the role of conduction and valence bands played by 
the upper and lower "Hubbard bands" with the "band 
gap" Amou- To make a pair with energy Amou, the dou- 
blon and hole must be separated from each other in a 
virtual process until they become real (on-shcll). The 
separation is on the order of A Mott /Fo an d the probabil- 
ity for this to happen is given by Green's function, i.e., 
V oc G(«AMott/^b) with k = ir/2. Because Green's func- 
tion decays exponentially in the Mott insulating phase 
[Eq. we obtain a production rate exponentially de- 

pendent on the electric field, i.e., Eq. ([3"5| . 

If we compare Eq. ([37| with Schwinger's threshold in 
QED [34|, HH, we notice that the correlation length £ 
plays the role of the Compton wavelength A = h/m e c. 
In the small U limit, correlation length (soliton length) 
is £ = 2t> e ff/AMott with v P a = 2 + U/2-n + ■ ■ ■ the speed 
of the charge excitations 1451 . In this limit, we recover 
the Landau- Zener result 1131] 
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(DC, Small U) (39) 



analogous to the Schwinger mechanism in QED with 
F th cx (gap) 2 Blll- 
In the large J7-limit, we have ~ \\\(U / gr) (g ~ 
4.3) [HI and production rate shows an interesting power 
law behavior 
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(DC, Large U) 



(40) 



with the hopping parameter r recovered. This result can 
be understood intuitively from a strong-coupling argu- 
ment. After pair creation in the Mott insulating ground 
state, the doublon must hop Amo« /Fo — U/ F sites away 
from the accompanying hole to become on-shell. The am- 
plitude decreases by a factor (■^■) for each hopping, and 
thus, we are led to Eq. (|4"0)) . In Fig. [5J we plot the U 
dependence of the threshold (Schwinger limit). 



IV. KELDYSH CROSSOVER IN AC FIELDS 

Next, we study a situation where a strong laser repre- 
sented by 
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FIG. 5: (color online) Schwinger limit (= tunneling threshold) 
[l5j of the ID Hubbard model at half-filling. In (c), the solid 
line is the Landau-Dykhne result given by Eq. (|35p . whereas 
the dashed line is its approximate form Eq. (|37[) . The Landau- 
Zener result in Eq. (|39l) with v e g = 2 is plotted as a dotted 
line. The Landau-Zener result is only accurate for a small U. 
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F(t) = F sin nt 



FIG. 6: (color online) Tunneling probability of the dh-pair 
obtained for the (7 = 8 Hubbard model in an AC fields, (a) 
and (c) correspond to -Fo = 1.0 and Q = 4.5, respectively, 
which are in the multi-photon regime, while (b) and (d) are 
for Fo = 1.0, and SI = 0.1, respectively, which are in the 
quantum tunneling regime. In (a) and (b), the tunneling 
probability is indicated by the size of the circle plotted on 
top of the dh-pair spectrum. 



(Fo: field strength, 57: photon energy) is applied to a 
Mott insulator. Experimentally, this models photocar- 
rier injection, which is the initial process in the photo 
induced insulator to metal phase transition [l], Q. In 
standard photocarrier injection, the laser's photon energy 
51 is set to the absorption peak, which is slightly above 
the Mott gap Amo«- However, herein, we are interested 
in the nonlinear process induced by subgap lasers, i.e., 
il < Amou- 

The tunneling probability V v as well as production rate 
T p can be calculated using Eq. (|2~4")) and Table HI We 
note that the sign in the Jacobian is determined so that 



(41) InxDpi^ > is satisfied. 
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FIG. 7: (color online) (a) Total production rate of the U — 8 Hubbard model in an AC field. The dashed line corresponds 
to Keldysh crossover 7 = f , where the excitation mechanism changes from multi-photon absorption to quantum tunneling. 
Schwinger limit is F t h = 1.668. (b), (c) Schematic pictures of the nonlinear excitation in the two regimes. In the quantum 
tunneling regime (Fig. (b)), the doublon-hole distribution becomes momentum independent, which means that the "upper 
Hubbard band" becomes populated by photocarriers. 



An interesting feature of photocarriers generated by 
nonlinear subgap lasers is that one can control the distri- 
bution of doublon-hole pairs by changing the photon en- 
ergy CI. In Fig. [51 we plot the momentum resolved tunnel- 
ing probability for the U = 8 Hubbard model. We notice 
that the distribution in the momentum space changes 
drastically when photon energy is changed. When Cl is 
large, the generated dh-pair is localized near the gap 
AE ~ Amou- The peak becomes broader as the field 
strength Fq becomes larger. On the other hand, when 
CI is small, the dh-pair becomes uniformly distributed 
in the p space. In the small CI limit, we approach the 



DC field case, where tunneling probability has no p- 
dependence, c. f., Eq. (f38| . In fact, the excitation mech- 
anisms in the two regimes are different. For small CI and 
large field strength, quantum tunneling is responsible for 
dh-pair creation. On the other hand, when CI is large, 
multi-photon absorption is the excitation mechanism. If 
we change photon energy and laser strength, there is a 
crossover between the two regimes, which is the Keldysh 
crossover [3l[ mentioned in the Introduction. We can di- 
rectly see this from the analytical expression of the p = 
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tunneling probability 



p=o ^ cxp 



2A Mo tt7 

y mi j 



/(7) 



(AC-fields), (42) 



7> 1, 



exp (-f^fffl -^ 7 2 + •••)) 7«1, 



(43) 



This result is obtained with the help of th e ap proximation 
in Eq. (|36[) (h ~ 1.47 and see footnote[52| for function 
/) . We note that this expression is identical to the QED 
result H?} with a redefinition of Keldysh's adiabaticity 
parameter (3l| 



(44) 



The crossover is characterized by the Keldysh line de- 
fined by 7 = 1. In the multi-photon absorption regime 
[7 <C 1; Fig. [TJc)], tunneling probability has a power law 
dependence on field strength. The power 2AMott/^ is 
twice the number of absorbed photons. As stated above, 
photocarriers are generated near the excitation gap. 

On the other hand, in the quantum tunneling regime 
[7> 1; Fig. [71(b)], tunneling probability shows a thresh- 
old behavior with an exponential suppression. In this 
regime, the photocarriers arc distributed almost equally 
in the momentum space. In Fig. [JJ (a) , we plot the total 
production rate as a function of photon energy and field 
strength. The Keldysh crossover line is indicated by a 
dashed line. In the quantum tunneling regime, the pro- 
duction rate quickly increases as field strength exceeds 
tunneling threshold, which is Fth = 1.668 for the U = 8 
Hubbard model. 

When we compare the present result with those of pre- 
vious studies on quantum tunneling in AC field back- 
grounds such as atom ionization [3l| and nonlinear 
QED [36W38} . we notice that the Keldysh crossover is 
quite universal and is not limited to the Hubbard model. 
Previous studies were conducted on non-interacting sys- 
tems where excitation occurs between single particle 
gaps, whereas in the present case, the system strongly 
interacts and the origin of the excitation gap is a many- 
body effect. The basic idea of the Keldysh crossover 
survives in many-body systems and expressions such as 
tunneling threshold (Schwinger limit) [Eq. ([3"T]) ] and the 
Keldysh's adiabaticity parameter (Eq. (|44j) ) are valid 
where the many-body features are renormalized on the 
correlation length £. However, if we carefully consider the 
long time dynamics, differences between non-interacting 
systems and many-body systems can be observed. This 
is examined below in the next subsection. 



A. Comparison with numerical results 

To examine the applicability of the Landau-Dykhne 
method, we compare its result with that of td-DMRG 
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FIG. 8: (color online) (a) Time evolution of doublon density 
in AC-fields calculated by td-DMRG for the U = 8 Hubbard 
model (L = 30, F = 4, Q, = 2). (b) Time evolution of the 
averaged doublon density d. Results for L = 30 and L — 20 
are plotted. Upper panel is the electric field F(t) = Fo cos fit. 
Doublon density shows intial increase and then saturation oc- 
curs, (c) Resonant oscillation between the ground state and 
the excited state with a neighboring doublon and hole pair. 



[511] . The ID Hubbard model on an L-sitc open chain 
with the Hamiltonian 



H(t) 



T + C jcr C j+l(T 



(45) 



3,<? 



is studied. This Hamiltonian is identical to the previ- 
ous one [Eq. ([5])] in the infinite size limit, and is related 
through a gauge transformation. We start time evolu- 
tion from the ground state |0) obtained by the finite-size 
method and apply the electric field for t > 0. In the cal- 
culation, the DMRG Hilbcrt space is m = 200, the time 
step is 0.01, and the length of the chain is L = 20, 30. 

There are several interesting features in the time evo- 
lution of the doublon density 



(46) 
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FIG. 9: (color online) Total production rate for the U = 8, 7 
Hubbard model with photon energy $7 = 1 plotted against 
field strength F ((b) is the logarithmic plot). The td-DMRG 
result is obtained using the fitting of Eq. ((48]), which is com- 
pared with the Landau-Dykhne result. Photon energy is 
Q, = 1 and system size is L = 30. 



as plotted in Fig. [8j First, in addition to the over- 
all increase, a density wave of doublons of the form 
di(t) ~ d(t) + (—l) l 5d(t) occurs. We note that the elec- 
tron density n.j = X) C r( n 'o') does not show such modu- 
lations. In addition, the finite-size effect is present in 
the calculation. This shows up in the doublon density 
as a difference between the edge and bulk values. To 
eliminate the finite-size effect as well as the density wave 
from the analysis, we define the averaged doublon den- 
sity d(t) as an average within the central region of the 
chain. For example, for an L = 30 system, the average 
is taken over the middle 10 sites. The time evolution of 
the averaged doublon density is plotted in Fig. [5Jb) for 
L = 20 and L = 30, and we see that size dependence is 
not large after averaging. The averaged doublon density 
shows a fast oscillation with a period of 2ir/U. One pos- 
sible explanation for the appearance of the density wave 
and time oscillation is the locality of dh-pair production. 
The correlation length £ gives the typical size of the dh- 
excitation created by electric fields. Because £ is short 
when U is large, e.g., £ < 1 when U > 10 [see Fig. EJb)], 
most of the dh-excitation takes place between neighbor- 
ing sites as in Fig. [5{c). The energy difference between 



the ground state this state is U, which leads to temporal 
oscillation. 

The averaged doublon density plotted in Fig. [8jb) 
shows an increase and the speed of increase becomes 
larger in stronger fields. After the initial increase, we 
notice that saturation takes place although the AC field 
is still present. If the field strength is sufficiently large, 
the long-term driven state is a metallic state called the 
"photo- induced Tomonaga-Luttinger-like liquid" [22| . In 
this state, it was numerically shown by calculating corre- 
lation functions that spin-charge separation takes place. 
However, in this study, we restrict ourselves to nonlin- 
ear dh-pair creation, which is responsible for the initial 
increase in doublon density. The production rate of the 
dh-pairs is obtained from numerical data by employing 
the fitting 

d(t) =d + a(tanh(fa) + 1) (47) 

and production rate is identified as the initial slope 

Tdmrg = ab, (48) 

In Fig[9] we plot the total production rate obtained by 
td-DMRG and compare it with the analytical result cal- 
culated from the Landau-Dykhne formula, i.e., Eq. (|29p . 
Although the fitting is difficult owing to the fast os- 
cillation of doublon density, the results seem to agree 
quite well. The numerical result agrees with the Landau- 
Dykhne result in the weak field regime up to Fq < 4. The 
deviation at large Fq is expected because the Landau- 
Dykhne formula is known to work only near the adiabatic 
limit, i.e., P<1. In the log- log plot, we clearly see the 
crossover from the weak field multi-photon absorption 
regime, with a power law behavior, to the quantum tun- 
neling regime. 

B. Single pulse 

Using the Landau-Dykhne method, we can extend the 
analysis to the case of a single pulse. There are two moti- 
vations. First, recent ultra fast pulse lasers are becoming 
so short that a single pulse excitation is now experimen- 
tally realizable. Second, we wish to further examine the 
validity of the Landau-Dykhne method from a theoretical 
perspective. In this regard, a single pulse is ideal because 
problems such as saturation and driven steady states do 
not occur, as is the case with a continuous AC field. The 
pulse we study here is 

F = F cosh -2 (i/cr) (49) 

and the Jacobian needed for the calculation is given in 
Table |TJ We set the duration to a — 2. 

In Fig. [TOTb). we plot the time evolution of doublon 
density for a ID U = 10 Hubbard model on an open 
L = 50 chain calculated by td-DMRG. Doublon den- 
sity initially increases and then oscillates with a period 
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FIG. 10: (color online) (a) Shape of pulse field with duration a — 2. (b) Time evolution of the doublon density d(t) for the 
U = 10 Hubbard model obtained by td-DMRG. (c) Increase in the doublon density Ad obtained by td-DMRG and the total 
production T obtained by the Landau-Dykhne method. See main text for the error bar. Inset: p-resolved distribution function 
of doublon-hole pairs obtained by the Landau-Dykhne method. 



roughly given by T = 2-k/U. This oscillation is similar to 
the AC field case seen in Fig. [8Kb) . Increase in the dou- 
blon density Ac? is plotted in Fig. ITUf c). where the error 
bars are determined by the oscillation width near t = 15. 
At Fo < 4, Ad shows a monotonic increase, whereas it 
shows an irregular decrease at Fq = 5. 

Let us compare the numerical result with Landau- 
Dykhne results. In the inset of Fig.fTOtc). the momentum 
resolved tunneling probability V v is plotted. The distri- 
bution has a peak at p = and broadens as the field 
becomes stronger. Eventually, as seen in the Fo = 3 
data, high-energy dh-pairs with momentum p ~ ±n be- 
come excited as well. We compare the total production 
rate T obtained from the Landau-Dykhne approach with 
Ad. Comparison with td-DMRG revealthat the Landau- 
Dykhne method gives the correct threshold behavior and 
is reliable even for relatively strong electric fields. How- 
ever, the Landau-Dykhne method does not capture the 
irregular drop at Fq = 5 and we detect a small deviation 
at Fo = 3. As mentioned above, the Landau-Dykhne 
method ignores the interference effect due to multiple 
quantum tunneling (4(|, which becomes important as pair 
annihilation of doublons and holes become activated. In 
fact, when the effect of interference is strong, quantum 
tunneling can be s upp ressed by dynamical localization in 
the energy space [23j|. Such an effect is expected to be- 
come important when the field is strong, and we think 
that it explains the differences between the numerical and 
Landau-Dykhne approach results. 



V. EXPERIMENTAL FEASIBILITY 

In this section, we discuss the experimental feasibil- 
ity of the nonlinear doublon excitations predicted in 



this study. Candidates of physical systems range from 
fcrmionic cold atoms to solid state crystals. For exam- 
ple, if we use fermionic cold atoms in an optical lattice, 
the present theory can be experimentally verified by in- 
ducing a time-dependent tilt of the lattice. This mimics 
the effect of an electric field and a direct measurement of 
the doublon density is possible as well [HI [H| • In solid 
states, we must estimate the threshold field strengths of 
candidate materials and compare them with the peak 
field strengths of present day lasers. In addition, because 
doublon density is not a directly observable quantity in 
solids, we must seek alternative detection methods. We 
discuss these issues below. 



Threshold of ID Mott insulators 



Pump-probe experiments with a terahertz (THz) laser 
are ideal setups to verify the Keldysh crossover described 
in section HV1 Current THz pulse lasers can be as strong 
as 1 MV/cm 0,0] and the typical photon energy is O = 
4 meV corresponding to 1 THz. 

On the material side, the candidate ID Mott insulator 
ranges from organic crystals to cuprates. The tunneling 
thresholds of several materials estimated by the Landau- 
Dykhne + Bethe ansatz method are shown in table [TTJ 
In the list, the material with the smallest threshold is 
ET-F 2 TCNQ (ET-salt). Let us discuss the possibility 
of generating dh-pairs in this material with a THz laser. 
The tunneling threshold (Schwinger limit) for the mate- 
rial is £"th = 3 MV/cm. When using a laser with peak 
strength -Ei ascr = 1 MV/cm @,0|, we are still below the 
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r(eV) U(eV) a(A) 


A M ott(eV) £(a) £th(MV/cm) 


ET-F 2 TCNQ 
[Ni(cnxn) 2 Br]Br 2 
Sr 2 Cu0 3 


0.1 1 10 
0.22 2.4 5 
0.52 3.1 4 


0.7 1.1 3 
1.6 1.0 16 
1.5 2.1 9 



TABLE II: Mott gap Ajviott, correlation length £, and tunneling threshold (Schwinger limit) _F t h = ea£" t h calculated by the 
Landau-Dykhne method with the Bethe ansatz. Material parameters (r: hopping, U : onsite repulsion, a: lattice constant) are 
obtained from Ref. @] (ET-F 2 TCN) and Ref. [H ([Ni(cnxn) 2 Br]Br 2 , Sr 2 Cu0 3 ). 



threshold, and the tunneling probability is 



V = CXp — 7T 



th 



-Ela 



X 10" 



(for ET-salt). 



(50) 



This is too small to trigger photo-induced metalliza- 
tion, and therefore, a stronger light source is needed. 
Quite recently, amplification of laser field strength us- 
ing a metamaterial structure has been proposed^, [55j . 
It is reported that peak strength can be as large as 
Smctamatorial =4 MV/cm || . If this technique can be 
applied, tunneling probability becomes 



exp 



Eth 



-^metamaterial 

0.1 (for ET-salt). 



(51) 



This means than one can perform a ~10% photodoping 
with a single pulse. This value is large enough that one 
can trigger a photo-induced phase transition. 



and the pump laser must exceed this strength. With 
stronger fields, more photocarriers are excited and the 
measurement is expected to be more feasible. 



C. Dielectric breakdown and nonlinear transport 



1 



metallic/ 



negative differential resistance 



threshold behavior 



- aV electric field F 



B. Angle resolved photoemission spectroscopy of 
the upper Hubbard band in the quantum tunneling 
regime 

As stated in Section IPVl [sec Fig. UJh), (c)], the dis- 
tribution of the produced dh-pairs changes drastically 
when the laser is shifted from the multi-photon absorp- 
tion regime to the quantum tunneling regime. In the 
quantum tunneling regime, high energy dh-pairs are gen- 
erated. This means that the entire "upper Hubbard 
band" becomes populated by photocarriers. Thus, using 
a strong THz laser as the pump, it is possible to study 
the band structure of the upper Hubbard band with a 
real time angle resolved photoemission spectroscopy tech- 
nique. A necessary condition is that the Keldysh param- 
eter (Eq. (gU) 



(52) 



is well below unity. From Table [TT1 the Keldysh crossover 
field strength for a 1 THz laser (0 = 4 meV) is 



E, 



crossover 



4 x 10" 2 MV/cm (for ET-salt), (53) 



FIG. 11: (color online) Typical /^/-characteristics of strongly 
correlated insulators. We have a threshold behavior, negative 
differential resistance and a transition to a metallic state. The 
present theory is only applicable to explain the threshold be- 
havior in the small current regime shown in the dashed box. 

In nonlinear transport experiments, a threshold be- 
havior in the TV-characteristics is found in many cor- 
related insulators. Materials range from Mott insulators 
[H, charge- ordered systems [H, and materials showing 
a neutral-ionic transition Q. In Fig Jill we plot typical 
JV-characteristics. In many cases, a threshold behavior 
as well as a region with negative differential resistance is 
present. Let us make a comment on the threshold behav- 
ior appearing in the small current regime, i.e., dielectric 
breakdown [region inside the dashed line in FigfTT]. 

Although doublon density is not a measurable quan- 
tity, we can relate the dh-pair production rate to current 
itself. Theoretical studies [UHil suggests that, in non- 
linear transport of Mott insulators, current has two major 
contributions 



J = aF + aT. 



(54) 



The first term is the standard linear response due to ther- 
mal carriers with a temperature dependence Othermai oc 
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e -A M ott/7\ The second term is proportional to the dh- 
pair production rate with a nonperturbative threshold 
form 



r = 2-exp(-^] 



(55) 



a is a nonunivcrsal proportionality constant which de- 
pends on the coupling to the electrode and other factors. 

One may think that photocarriers induced by dh-pair 
production may contribute to linear response, i.e., the 
term oFq. This is not true in dielectric breakdown oc- 
curing in the quantum tunneling regime. The reason 
is because the dh-pairs are in an "infinite temperature 
state" [56| . As shown in Section IIII1 the distribution 
of dh-pairs has no momentum-, and thus energy-, de- 
pendences. This corresponds to an infinite temperature 
state, i.e., e ~ E / kBT with T — > 00. The conductivity of 
an infinite temperature state is zero, and therefore, there 
are no linear rcponse contribution in the quantum tun- 
neling regime. Instead, the current is dominated by the 
second term in Eq. (|54[) which is proportional to the dh- 
pair production rate. The doublons and holes are pair 
created, separated from each other by the electric field, 
and are measured as current when they are absorbed by 
electrodes. This feature is consistent with numerical re- 
sults obtained in a static system coupled to electrodes 

After dielectric breakdown, the /^-characteristics 
show interesting nonlinear behaviors. Although this is 
far beyond the applicability of the present theory, let us 
consider existing literatures. If voltage drop is measured 
as a function of current, there is a regime where nega- 
tive differential resistance is realized [J-Q. The origin 
of negative differential resistance is not fully understood 
yet. It was pointed out, with a careful comparison with 
experimental data, that the temperature increase of the 
sample due to Joule heating can explain it Q. They 
suggest that negative differential resistance occurs when 
the temperature dependence of conductivity is large. A 
related theoretical paper explained negative differential 
resistance in disordered films via the heating mechanism 
[54| . A more dramatic proposal is based on nonequi- 
librium first-order phase transition proposed by Ajisaka 
et. al. [25|, where negative differential resistance is ex- 
plained through a phase bi-stability. Negative differen- 
tial resistance was also found in a model in high energy 
physics, namely the supersymmetric QCD in the large 
iV" limit J2(|. A microscopic understanding of the non- 
linear transport properties of correlated systems from a 
universal viewpoint is an interesting callenge. 



D. Optical sum rule 

From Eq. (|54j) , we can derive an interesting relation- 
ship between the optical sum rule and the doublon pro- 
duction rate in DC-electric fields. Here, we consider the 
low temperature case where contributions from thermal 



carriers are negligible, and the field strength is below the 
threshold [region inside dashed line in Figure [TT] . 

The optical sum rule is commonly used by experimen- 
talists as a means to "measure" carrier density from the 
absorption spectrum[l|. The optical sum rule for the 
Hubbard model (e.g., 49|) states that the frequency in- 
tegral of the absorption spectrum, which we call N e g fol- 
lowing Ref. [l|, is related to kinetic energy as follows: 



dw 



o x (w) 



2L 



-K, 



(56) 



Here, K is the time average of the expectation value of 
the kinetic term in the Hamiltonian, i.e., 



K(t) = (*(t)| - r^2(e^c] +la c ja + c] a c 3+la )\^{t)) . 



(57) 



Below the threshold, dh-pair creation is a rare event and 
current is very small. In this regime, energy dissipation 
by external degrees of freedom, e.g., phonons, is negligi- 
ble and we can use Joule's relation 



j t (H{t)) = LJF 
to relate change in energy 

(H(t)) = K(t) + ULd(t) 



(58) 



(59) 



to the current J. Using Eq. (|54|) as well as 4id(t) = T, 
we obtain 



dt 



N ce = (U- aF )T. 



(60) 



This formula states that, in static electric fields, optical 
sum N e g increases linearly in time and the speed is pro- 
portional to the dh-pair production rate. This relation- 
ship opens a way to measure doubon increase by optical 
experiments. 



VI. 



CONCLUSION 



We developed an analytical theory for nonlinear pair 
excitations of doublons and holes in a ID Mott insulator 
subject to DC, AC and pulse electric fields. The theory 
is based on the Landau-Dykhne method combined with 
the Bethc ansatz. In an AC-field, the theory predicts a 
crossover between multi-photon absorption and quantum 
tunneling when the strength and photon energy of the 
field changes. Comparison with numerical results by td- 
DMRG shows that the analytical theory is reliable up to 
moderate field strength. 

There are several limitations in our theory. Perhaps 
one of the most important open issues is the treatment of 
temperature effects. It is unclear if the tunneling prob- 
ability has a direct temperature dependence. Numeri- 
cal results obtained by nonequilibrium DMFT suggest 
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no or very small temperature dependence in the quan- 
tum tunneling regime [l7|, whereas strong temperature 
dependence of the threshold is seen in a dielectric break- 
down experiment Q. Impact ionization, an avalanche 
like cascade growth of carriers due to field induced accel- 
eration, may be important in understanding these exper- 
iments. The origin of negative differential resistance and 
the properties of the nonequilibrium steady state is an- 
other important open problem ( Section [V C|) . We think 
that trying to answer these problems will lead to impor- 
tant innovations in nonequilibrium manybody physics. 

We acknowledge Philipp Werner, Kunio Ishida, David 
Pekker, Rajdeep Sensarma, Li Gao, Stuart Parkin, Ger- 
ald Dunne and Eugene Demler for valuable discussions. 
TO acknowledges support from Grant-in- Aid for Young 
Scientists (B), CUA and ITAMP. 



Note added: After the submission of the initial ver- 
sion of the manuscript, Lenarcic and Prclovsck published 
an interesting paper [531 ■ They studied the dielectric 
breakdown in a spin polarized Mott insulator and found 
that the threshold has a F t h oc A 3 / 2 dependence, which 
is different from the conventional Landau-Zener form. 
They pointed out that the origin of this difference comes 
from excitation dispersion. In the half-filled case, the 
dispersion is rclativistic cok oc 

whereas it is parabolic uik oc fc 2 + k 2 in the spin polarized 
model. Using Eq. ([55)1 . we can recover the F t ^ oc A 3 / 2 
behavior in the spin polarized case. In another article, 
a variant of Eq. (f35j) was used to study excitations in 
the attractive Hubbard model (5^- These examples show 
the wide applicability of the Landau-Dykhne approach 
in many-body problems. 
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